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Abstract 

Using the connection between intersection theory on the Deligne- 
Mumford spaces M 9jn and the edge scaling of the GUE matrix model 
(see fl2| , |l4}|), we express the n-point functions for the intersection 
numbers as n-dimensional error-function-type integrals and also give 
a derivation of Witten's KdV equations using the higher Fay identities 
of Adler, Shiota, and van Moerbeke. 



1 Introduction 

1 . 1 Overview 
1.1.1 

This paper is a continuation of |L2| and is very closely related to |14|. It 



was observed in [I2| and conceptually explained in [TJJ that the intersection 
theory on the moduli space of curves is closely connected to the edge scaling 
of the standard GUE matrix model (which can be called the edge-of-the- 
spectrum matrix model). Leaving the detailed discussion of this phenomenon 
to |14]], we obtain here some formal consequences of this connection, most 
importantly, an error-function-type integral formula for certain generation 
functions for the intersection numbers known as the n-point functions. We 
also show how to derive, using results of Adler, Shiota, and van Moerbeke 
UK, the KdV equations from this matrix model. 
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1.1.2 



Let M 9i „ be the stable compactification of the moduli space of genus g curves 
with n marked points. Let ipi be the first Chern class of the line bundle whose 
fiber over each pointed stable curve is the cotangent line at the zth point. 
We use the standard notation 

(r dl ...r dn )= [_ y)di = 3<7-3 + n, (1.1) 

<>M B ,n 

for the intersection numbers of these ^-classes. 

The celebrated conjecture of Witten |17| says that for the following gen- 
erating function 

fmi,...)= e <^--->n|' 

(fco,fci,...) 

the exponential exp(F) is a r-function for the KdV hierarchy in variables 

U 



T, 



2i+l 



(2i + l)!! 



This conjecture was inspired by matrix models of 2-dimensional quantum 
gravity, see for example j|, [5] for a survey. 



The proof of Witten's conjecture given by Kontsevich [TT[] uses another 
matrix model interpretation of the generating function F. 

1.1.3 

The purpose of this paper is to exploit yet another random matrix connection 



T2| to evaluate a different generating function for ( |1 . 1| ) , called the n-point 
function, in a closed form. 

Definition 1.1. We call the following generating function 

oo 

9^(xi, • • • , X n ^j ^ '^gi^'y. • • • , X n ) , 
9=0 

where 

? g (xi, ...,X n )= E ( T di ■ ■ ■ T dn) II X ^ ' 

E di=3g-3+n 

the n-point function. 
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The 2 and 3-point function were computed by Dijgraaf and Zagier, re- 
spectively. Several specialization of the 3-point function can be found in the 
paper ||. Our main result, see Theorem [I] below, is that the n-point function 
J is a certain specific multivariate error function (and hence, in particular, 
a multivariate hypergeometric function in the sense of Gelfand, Kapranov, 
and Zelevinsky, see Section 1.2.3| ). 



1.1.4 

We recall, referring the reader to 0, [14| for details, that the n-point functions 



turn out to be connected to the asymptotics of the following averages over 
N x N Hermitian matrices 



oo 



;i.2) 



\k=l 



N 



with respect to the standard Gaussian measure. Here [ • ] stands for the in- 
teger part. The asymptotics of ( |1.2|) is equivalent to the knowledge of the 
distribution of eigenvalues of a random matrix near the edge of the Wigner 
semicircle. This distribution is well known to be described by the Airy en- 
semble, see for example JT5 |. We also note that this Airy behavior near the 



edge of the spectrum is very universal for random matrix ensembles, see for 
example fL5|. 



1.1.5 

By a standard application of the Wick formula, the asymptotics ( |1.2j ) is 
equivalent to the asymptotics in the following combinatorial enumeration 
problem. Consider a surface E and a map on it with n cells, that is, a way 
to glue E out of n polygons by identifying the edges of the polygons in pairs. 
The asymptotics of the number of ways this can be done, as n is fixed and the 
perimeters of the polygons go to infinity at certain relative rates, is encoded 
in the asymptotics of (\1.2\) . A typical map with 3 cells on a genus g — 2 
surface E is shown in Figure |l|. Here one portion of the map is magnified 
so that to show the dendriform pattern that the map forms because the 
polygons are allowed to be glued to itself. In fact, the overwhelming part of 
the perimeter is typically contained in these trees as the perimeter goes to 
infinity. The macroscopic data of the map in Figure [I] is a trivalent graph 
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Figure 1: A typical map of large perimeter 



embedded in £ which is the same data as in HIT 



sec 



12L 111 for details of 



the connection of this approach to [11 



The asymptotics in this map enumeration problem is described by a cer- 
tain function map s (:ri, . . . ,x n ), see Section 2.1.7 of 
tions and Section |2.1| below for a formula for map s . 
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for precise defini- 
When the surface S 
is connected of genus g, we write map 9 instead of map s . Obviously, the 
asymptotics map s is multiplicative in connected components. 

There exist a number of other combinatorial asymptotics equivalent to 
the asymptotics (|1.2f ), such as the distribution of increasing subsequences in 



a random permutation, see for example 



1.1.6 

Another product of our analysis is an alternative derivation of Witten's KdV 
equations using the edge matrix model. It is based on the Fay identities 
methods developed by Adler, Shiota, and van Moerbeke in [|TJ. 



1.2 Formula for the n-point function 
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1.2.1 The function £ 

The key ingredient in our formula for the n-point function will be the follow- 
ing function of X\, . . . , x n 



£(xi, . . . , x n ) 



1 exp x fj 



2 n 7T n/2 Y\ y/X~ 



X 



/ ds exp ( — 

J Si >o V i=l 



4:Xi 



Si + Si+i 



;i.s) 



where the integral is over M> and s n+ i = Si. 

This integral admits a nice probabilistic interpretation, namely 

where the integral is over the space of nonnegative piecewise linear function 
of the form shown in Figure subject to the bridge condition s n+1 = s\. 




Xj X 2 



JC ? 



Figure 2: A piecewise linear function 7 



The measure W is the natural Gaussian measure: the increments of 7 are 
independent normal variables with variance equal to twice the length of the 
interval. 
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1.2.2 Main result 

The function £ is clearly invariant under a cyclic shift of the x^s. We want 
to have something symmetric in the X{S, so we set 

e°(x) = £(x ct( i),...,x (T ( S )), 

o-GS(n)/(12...n) 

where the summation is over coset representatives modulo the cyclic group 
generated by the permutation (12 . . . n). 

Let il n denote the set of all partitions a of the set {1, . . . , n} into disjoint 
union of subsets. For any partition a G Il„, with I = £(a) blocks, let x a 
denote the vector of size I formed by sums of Xi over the blocks of a. For 
example, if 

a = {l,3}U{2} en 3 

then £(a) = 2 and x a = (x\ + 23, X2). 
By definition, we set 

9(x u ...,x n )= J](-l)^ +1 £ (x a ). 

aen„ 

For example, 

9{x 1 ,x 2 ) = £(xi + x 2 ) - E(xi,x 2 ) , 

and, similarly, 

9(xi,x 2 , x 3 ) = £(xi + x 2 + x 3 ) - £(xi + x 2 , x 3 ) - £(xi + x 3 , x 2 ) 

- £(x 2 + x 3 , x x ) + E(xt, x 2 , x 3 ) + £(xi, x 3 , x 2 ) . 

In this notation, our result is 

Theorem 1. 




We remark that the right-hand side of ( |1.4|) makes sense for all g > and 
n > 1, even for the pairs 

G7,n) = (0,l),(0,2) 

for which the corresponding moduli space is problematic. The corresponding 
terms, however, are not polynomial in x, see Section [2. 6. 3 . 
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1.2.3 

The integral (|1.3|) fits inside a more general class of integrals of the form 



HQ) 



[ e Q(s) ds, degQ = 2. 



Such a multivariate analog of the error function function of coefficients 

of the polynomial Q, a multivariate hypergeometric function in the sense of 
Gelfand, Kapranov, and Zelevinsky. In fact, this is true for any polynomial 



Q 



E 



m=(rai,...,m„) 

because the integral I(Q) satisfies the obvious equations 

(n^-nsy '«>=». E^E* W 

and integration by parts gives 

I(Q) = ~ [ e Q ^ Si ^-Qds, 

which is equivalent to the homogeneity equations 



{y^rma m ^- + l\ I(Q) =0. 



for i = 1, . . . , n. 



2 Proof of the n-point function formula 



2.1 



We will assume that the reader is familiar with Section 2 of [I2"|. The main 



object we need from |12| is the function map 9 (x) which describes the N — > oo 
asymptotics of the number of genus g maps with n cells with perimeters 
~ Nx%, • • • , ~ Nx n , respectively. 
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We have the following formula which is Theorem 3 in fT^ 



e ^ map (af 



dx 



x 



2le(r)|/2 



— y 1 n 

)/2-l I An+mi 11 



(2.1) 



where the summation is over all trivalent ribbon graphs T of genus g with 
n cells, the product is over all edges e of T, z\ >e and z 2>e are the Zj's corre- 
sponding to the two sides of edge e, and 

|e(r)| =6#-6 + 3n 



is the number of edges of any graph T G T 3 



g,n- 



On the other hand, we have the following formula of Kontsevich which is 



the unique boxed formula in |l I 
1 



n 



r 6 r»„ ' Aut < r >l «") Zl -' + z ^ 



2 -H«(r>l/3 £ (r A ...r*,>n 

J] di=3g-3+n 



(2di 



Zi 2 *- 1 . (2.2) 



8=1 



We refer the reader to \TA\ for a very detailed discussion and proof of this 
formula. 



From ( p.lD and (|2.2j ), using the formula 



( " ' (II i f/// • -] S - m -- 



(2m)! 

7T — — S 

2 2m m! 



and the fact that the function map 9 is homogeneous of degree 3g — 3 + 3n/2, 
we obtain the following: 



Proposition 2.1. 

g \P^1 ) ■ ■ ■ j •fcr, 



7r n / 2 map (2xi, . . . , 2x n ) 



29 



11* 



1/2 



(2.3) 
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2.2 

The function map„ can be obtained as the Laplace transform of a certain 



natural object. Introduce the Airy kernel, see for example |L6 

Ai(ar) M\y)-M'{x)M{y) 



K(x,y) 



x-y 



where Ai(x) is the classical Airy function. 

The kernel K differs by scaling of variables from the kernel K used in |12| . 
It is natural to use K when working with map g (2x) instead of map 5 (x). We 
now introduce the corresponding modifications of the functions from Section 



2.1.8 of 12 



Let the function R be the following Laplace transform 

R(6,--.,0= / det[K(x t , Xj ))]dx. (2.4) 



By definition, set 



H(zi,...,x n ) = R «>> ( 2 - 5 ) 
aen n 



where Il n denotes the set of all partitions of the set {1, . . . ,n} into disjoint 
union of subsets, and x a is the vector formed by sums of Xj over blocks of a. 
For example, 

H(xi, x 2 , x 3 ) = R(xi, x 2 , x 3 ) + R(xx + x 2 , x 3 ) + 

R(xi + x 3 , x 2 ) + R(x 2 + x 3 , xx) + R(x l +x 2 + x 3 ) . 



Finally, set 



G(x u ...,x n )= ^ H(xj) ie5 H(xi)^ 5 , (2.6) 

5c{l,...,n} 



where the summation is over all subsets S and H(xi)nzs denotes the function 
H in variables Xj, i G S. For example 

G(xi, x 2 ) — 2H(x l5 x 2 ) + 2H(xi)H(x 2 ) . 

The function map 9 has a natural extension to disconnected surfaces S 
by multiplicativity. We have the following formula, see the last formula in 
Section 2.1.8 of @, 
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Proposition 2.2. We have 

G(xi, . . . , x n ) = ^ map E (2xi, . . . , 2x n ) , (2.7) 

where the summation is over all orientable surfaces E ; including disconnected 
ones. 

For example, 

G(xi,x 2 ) = 2jmap s (2x 1 , 2x 2 ) + map qi (2xi) map^ (2x 2 ) ■ 

9 91 ,92 

2.3 

Definition 2.3. Let Q(xi, . . . , x n ), n — 1, 2, . . . , be a sequence of function 
which are symmetric in their arguments. We define 

Q°(x 1 ) = Q(x 1 ) 1 
Q°(x 1 ,x 2 ) = Q{xi,x 2 ) - Q{x 1 )Q{x 2 ) , 
Q°(xi, x 2 , x 3 ) = Q(xi,x 2 , x 3 ) - Q{xi)Q(x 2 , x 3 ) - Q(x 2 )Q{xi, x 3 ) 
- Q(x 3 )Q(xi,x 2 ) + 2Q(xi)Q(x 2 )Q(x 3 ) , 

and so on, namely, terms of degree d in Q come with the coefficient (— l) d ~ 1 (d— 
1)! which is the well known Mobius function for the partially ordered set Il n 
of partitions of {1, . . . ,n}. We call these functions Q° the connected part of 
Q. 

It is clear that the connected part of the sum (|2.7|) is the following 
G°(xi, . . . , x n ) = map g (2xi, . . . , 2x n ) . 

9 

It is clear from ( |2.3| ) that we have 
Proposition 2.4. 

2 n/2-l n n/2 

3(xu ■ ■ ■ . x n ) = 1/2 G 

U x i 

Now it only remains to check that 



(2V3) ' 

(2.8) 
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2.4 



The sum G can be interpreted in a diagrammatic way as follows. We take 
our variables {x{\ and divide them first in two large groups in all possible 
ways as in ( |2.6| ). Next we divide each group into smaller subgroups, which 
are the blocks of the partition a in fl2.5|) . These blocks are then grouped 
together into the cycles of permutations which appear in the determinant in 

flU). 



The contribution of each such diagram to G is a product over the con- 
nected pieces of the diagram. Therefore, in the function G° every discon- 
nected diagram will cancel out and G° will be the sum over connected dia- 
grams only. 

The corresponding summands are the following. First, in (|2.6| ) every 
nontrivial subset S leads to a disconnected diagram and so the connected 
part of G is just twice the connected part of H. After that, in (|2.5| ) we can 



take any partition a, but then in ( |2.4| ) we need our permutation to have only 
one cycle. 

Denote the contribution of one long cycle to (|2.4|) by 



„ n 



(zi, z i+1 ) dz , (2.9) 



with the understanding that x n+ \ = x\. Note the invariance of the function 
(|2.9j) under the cyclic permutation of variables. 



The connected part of R is therefore 

R (x) = (-1)"' +1 E(x a(1) ,...,x CT(s) ), 

aeS{n)/(12...n) 

where the summation is over representatives of the right cosets of the sym- 
metric group modulo the cyclic group generated by the permutation (12 ... n). 
With this notation, we have the following 



Proposition 2.5. 



«en„ 



For example, 

G°(a;i, x 2 ) = 2R°(x 1 , x 2 ) + 2R°( Xl + x 2 ) = 2E( Xl + x 2 ) - 2E(x l5 x 2 ) 
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Similarly 



G°(x 1 ,x 2 , x 3 ) = 2E(xi + x 2 + x 3 ) - 2E(xi + x 2 , x 3 ) - 2E(x 1 + x 3 , x 2 ) 

- 2E(x 2 + x 3 , xi) + 2E(xi, x 2 , x 3 ) + 2E(xi, x 3 , x 2 ) . 

2.5 

Now it only remains to prove that 

£ = E. 

Consider the function ( |2.9| ) more closely. We have the following formula, 



see formula (4.5) in [jig] , 

poo 

K(z,w)= / Ai(z + a) Ai(w + a) da . (2.10) 
Jo 



Substituting this into ( |2.9| ) and interchanging the order of integration we 
obtain the integrals of the form considered in the following 

Lemma 2.6. 

f° xz w \ a • / U s , 1 fx 3 a + b (a-b) 2 \ 

L 6 + a) Al( ^ + h) dZ = 2^=x 6XP ( 12 - — X - 4T" ) 



Proof. Denote by / the function / = Ai(z + a) Ai(z + 6). The Airy function 
is a solution of the Airy differential equation Ai"(2:) = zAi(z). It follows that 

fzzzz - (42 + 2a + 26)/ M - 6/ 2 + (a - 6) 2 / = . 

This translates into a first order ODE for g = J e xz f dz. 

fx 2 _a + b 1 {a-b) 2 \ 
9x ~ \ A~^2T~2x~ + Ax 2 ) 9, 

from which it follows that 

h(a,b) fx 3 a + b (a - b) 2 \ 
9 = ^T eXP [-2-— X -^x-) ' 

where h(a, b) is some function of a and 6. The obvious equations 
faa-(z + a)f = Q, f bb -(z + b)f = 



12 



translate into the equations 

9aa -ag = g bb - bg = g x , 

which imply that h is constant. The fact that h = — = follows by taking the 
Laplace transform of ( [2.10D with z = w and comparing it with the known 



formula for it, see Section 2.6 of [|T2|. □ 
This lemma concludes the proof of the theorem. 

2.6 Examples and remarks 
2.6.1 

In the simplest example n = 1, there is no quadratic term in the exponential 
and we obtain the formula 

which corresponds to the well-known result 

. 1 rv , . exp(a; 3 /24) 

24 9 g\ x z 

2.6.2 

The case n = 2 reduces to the ordinary error function 

2 f x _ 2 
erf(x) = —j= I e u du , 

V 71 " Jo 

as follows. The following integral 



exp(— P(a + b) — Q{a — b) 2 ) da db 



o 



^ eX p (^-) (l - erf f ) ) ( 2.11 ) 



2Pv^ ^ V 

can be computed by introducing the variables u = a + b and v = a — b and 
integrating first in u and then in v. 
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The integral in ( |1.3| ) for n = 2 becomes the integral ( ^.11| ) with the 
following choice of parameters 

p _ Xi + X 2 q _ Xi + x 2 

2 " 4xiX 2 

Therefore 

3 



1 exp ((xi + x 2 y/12 



£(xi, x 2 ) = -, ■ 1 - erf -^x l x 2 {x 1 + x 2 \ 

2^/ix (xi+x 2 y/ 2 \ \2 

It follows that 

S(xi, x 2 ) = £(xi + x 2 ) - £(xi, x 2 ) 

1 exp((x 1 + x 2 ) 3 /12) (\ 



2v^F (xi+a; 2 ) 3 / 2 V2 
Now we use the expansion 



erf [ -\J xix 2 (xi + x 2 ) . (2.12) 



;* 2 / 4 <Tf(.r 2) ^=> 7—7— — tt x 2k+] 



k=0 



which can be proved, for example, by noting that the function erf(x/2) satis- 
fies the equation u xx + 1 u x = and hence the function e x / 4 erf (x/2) satisfies 
the equation u xx — \ u x — = 0. It follows that 



I / x 3 x 3\ °^ fci (1 \ k 

^ X2) = ^Tx~ 2 6XP U + to (2*Tl)! U XlX2(Xl + X2) J 

(2.13) 

By the string equation 

^(xi, x 2 ,0) 

J(X!, £ 2 ) = ; , 

Xi + x 2 

and so again we find agreement with the known formula for J(xx, x 2 , 0), see 
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2.6.3 

Observe that the only non-polynomial term in ( [2.13 ), which is comes 



Xi + x 2 

from the exceptional case (g,n) = (0,2). 

Also observe that the term £(xi + x 2 ) cancels out completely in ( |2.12| ) 
with a part of £(xi, x 2 ). This cancellation can be seen a-priory as follows: all 
coefficient of "J are rational numbers and the term £(x\ + x 2 ) has the wrong 
power of 7r. This argument works in general to identify terms making no 
contribution to S- 

In fact, the formula (|1.4j ) is full of cancellations of this type. Several 
criteria for identifying irrelevant terms will be developed in the next section 
where we will be dealing with the KdV equations. 



3 KdV hierarchy 

3.1 Strategy 
3.1.1 

In this section we will give a derivation of the Witten's KdV equations which 
uses our current machinery and the Fay identities techniques developed by 
Adler, Shiota, and van Moerbeke in |]J . Namely, we will prove that the series 
expansion of &(x) about x = is given by coefficients of one specific KdV 
r-function. See []11[ , and also for example M, for the exposition of how this 
was established originally using Kontsevich's matrix model. 

3.1.2 

We will need some qualitative facts about the series expansion of S(x) about 
x = 0, the first one being that all exponents in that expansion are half- 
integers. Consequently, after sufficiently many differentiations, every terms 
in this expansion blows up as Xj — > +0 for any i = 1, . . . , n. In other words, 
the coefficients of S(x) can be determined from singularities of 9 and its 
derivatives as — > +0. 
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3.1.3 



Proposition |2.5| and (2.8) express the function 9 in in terms of the function 



R which, by definition, is the following Laplace transform 

R(6,...,60= / e<-Mdet[K(xi, Xj )]dz. (3.1) 

Hence the x — > singularities of 9 and its derivatives are determined by the 
similar singularities of R. Clearly, the singularities of R and its derivatives 
as & — > +0, are determined by the x — > — oo asymptotics of the integrand in 
dp). 

The x — > — oo asymptotics of the integrand in (|3. 1|) can be computed 
from the classical asymptotics of the Airy function, however it is rather com- 
plicated and difficult to use in ( |3.1| ) directly. Fortunately, there are several 
qualitative criteria, a sort of selection rules, which help identify many terms 
in the asymptotics of det [K( irrelevant. 

3.1.4 

Observe that because all exponents in the expansion of 9(#) are half-integers, 
no information is lost by computing — ^-j 9(%) instead of 9- Indeed, 
the operator 

i-^e End c (VxT^Clxf, • • • , xf]) 

has no kernel. On the other side of the Laplace transform ( |3.l| ), this means 
that we are free to multiply the integrand by any factor of the form (xi —Xj). 

This is convenient for the following reason. The expression of 9 in terms 
of R involves terms which are Laplace transforms of distributions supported 
on diagonals like the E(x\ + x-z) term in (p.!2[) . Also, the asymptotics of 



the integrand in (R.ll) contains denominators of the form which require 

special handling resulting in certain contributions to the asymptotics from 
the diagonals. Since these contributions disappear after multiplying by ap- 
propriate factors of the form (xi — Xj), they always balance out in the final 
answer like they did in ( |2.12| ). 

In other words, this argument shows that if a term in the asymptotics of 
det [K(xj, X,)] becomes negligible after multiplication by (x^ — Xj), then it 
can be discarded. 
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We will call such negligible terms the diagonal terms. For example, The- 
orem |l| can be restated as saying that 



9(xi, . . . , x n ) = —2 £° \ 7^) + diagonal terms . 

n x i 

Another qualitative way to eliminate the diagonal contributions is to observe 



that they come with the wrong power of ir, see Section 2.6.3. 



3.1.5 

Another category of negligible terms are the oscillating terms in the asymp- 
totics of det [K(xj, Xj)] . The x — * — oo asymptotics of the Airy function 
Ai(x) is oscillating and most of the terms in the asymptotics of the integrand 

in ( |3.1|) will have factors like exp (j^ix^/j for some k = 1, . . . ,n. We now 
observe that such a term will never blow up as — > +0, even if we multiply 
it by an arbitrary large power of Xk- Indeed, we can assume that we already 



got rid of all denominators as explained in |3.1.4| above, and then it is enough 
to show that 



e -^+» 3/ V dx = o{i), £ ^ +o 



where c > and a is arbitrary. After a change of variables, it becomes 
equivalent to 

e -t*>/* e ix x a dx = ; ^ _> +0 , 

for some other constants a and c. Now, if we integrate by parts integrating 
e tx and differentiating the rest, we can decrease a so that to make the integral 
with £ = absolutely converging, thus proving our assertion. 

3.1.6 

Finally, we are interested in terms in S(x) of positive degree in all Xi. Recall 
that all terms of 5(x) have positive degree in all Xj's with the exception of 
the g = terms for n = 1,2. 

3.2 Asymptotics of the Airy kernel 
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3.2.1 

Introduce the following functions 



a[x) = > cifcX 

fc=0 



-3k 



a k x 



-Sk+l 



k=0 



where a = a = 1 and 



vfc (6/s-l)!! . l + 6k 
-i) , a k = -afc. 



72 k (2k)\ ' ~ re 1 - 6A; 
Also set, by definition, 

A(x) = exp (^f^ 3 ) , A(x) = exp (ifa; 3 ) a(x) . 

and 

q(x) dfa) - a(y) a(x) A(x) A(y) - A{y) A{x) 

a(x,y) = , A{x,y)- 



x — y x — y 

We have (see e.g. @) the following x — > +oo asymptotics 

1 



Ai(-x) 
Ai'(-x) 



2 x /vFix 1 /4 
1 

~2v / 7ra; 1 /4 



/ 4 A (x 1 / 2 ) - e^M (-x 1 / 2 ) 
/ 4 i (x 1 / 2 ) - e~^ 4 i (-x 1 / 2 ) 



This asymptotics remains valid for complex x such that | argx| < |7r. 
3.2.2 

It follows that away from the diagonal x = y we have 



y e,e'=±l 



ey/x + e'^fy 



where the sum is over 4 possible combinations of the signs e and e'. Because 
the function K(— x, —y) is analytic at x = y and its asymptotics remains valid 
in the complex domain, this asymptotics can be extended, using contour 
integrals, to the diagonal x = y. 
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By multinearity of the determinant, we have 



det [K(-x i ,-y j )] 1 



TJTTJl det 



e,e'G{±l} n 



A (si-s/Xi, Ejy/T/jj 
X% Uj 



(3.2) 



where the summation is over 4 n choices of signs vectors e and e' . Observe 
that by definition of A we have 



det 



A(si,3/j) 
Xi + yj 



exp (fiE^ 3 + V?)) det 



a (xj,yj) 
Xi + yj 



This last determinant can be simplified using the higher Fay identities |l 
which will be discussed in the following section. 



3.3 Fay identities and r-functions 
3.3.1 

In this section we collect, for the reader's convenience, some background 
material about T-functions and Fay identities QXQ . We will slightly deviate 
from the notational conventions of the book || by V. Kac. For our purposes, 
it will be convenient to use the notations of the Appendix to |13|] . 



3.3.2 

Given a matrix M 6 GL(oo), we will denote by (M) the charge vacuum 
matrix element 

(M) = (Mv h v % ) 

for the action of GL(oo) in the infinite wedge space. For any M, the function 

r M (t) = (T + (t)M) , 

of the variables t = (t\, £2, • • • ), is a r-function for the KP hierarchy. In our 
case, we will additionally assume that M is upper triangular, which implies 
that M*v = t>0 and hence 

tm(0) = 1 . 
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3.3.3 

Consider the following matrix element 



vi=l 



where ip{z) and ip*(w) are the standard generating functions for the fermionic 
operators. On the one hand we have (see Corollary 14.10 in || reproduced 
in (A. 14) in |3) 



i/j(z) ip*{w) 



'zw 



z — w 



F-({z} - {w}) r+dw- 1 } - {z- 1 }) 



where 



z^ z 

{z}=[z,-,j,... ) . 



From this and the commutation rule 

r + (t)r_( s ) = e SfcifcSfe r_( s ) r+(t) 

it follows that 

A(z)A(-w) 



T 



(EK 1 } - E{% 1 » 



(3.3) 



no* - 

where A denotes the Vandermonde determinant and Ho = Yli y/zjuk- 
3.3.4 

On the other hand, consider the operators iPm{z) = M~~ 1 ip(z)M, and simi- 
larly ip* M {w). Using the commutation rules and M*v% = v$, we compute 



det [ (ipM(zi) ip* M {wj)) ] = n det 



r {{w; 1 } - {zj 1 }) 



Zi - Wj 



• (3-4) 



Here the second equality is based on Wick's theorem, or equivalently, fol- 
lows from the observation that, acting on the vacuum, the operators ipM^Wj) 
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remove vectors k, where k G {— |, — |, . . . }, and then the operators ipM{zi) 
have to put all these removed vectors back, in all possible orders. 

Combining ( |3.3| ) with ( |3.4|) and interchanging the roles of z and w, we 
obtain the identity 



det 



r a* 1 } - {wj 1 }) 



Zi - Wa 



A(z)A(-w) 



(EKVEK 1 }), (3.5) 



which is the form of the higher Fay identities p| that we will need here. 



3.3.5 

From now on, we will be interested in one particular r-function, see e.g. 
corresponding to the matrix Ma which acts as follows 



m 



Ma k 



^2 di k + 3i . 
eg k + 3i . 



k + 1/2 is even . 
k+1/2 is odd. 



That is, the matrix M looks like this: 



M A 



a 



a i 



a 



ai 



a 



«2 



Oi 



a 



Ct>2 



a i 



where the empty spaces represent zeros. This matrix commutes with the 
bosonic operators otu for k even and hence the r-function 



ta = t Ma 



does not depend on tk with k even, meaning that it is a r-function for the 
KdV hierarchy. It also follows that 



r A {t - {w}) = r A {t + {-w}) 



(3.6) 



for any w. 
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3.3.6 

It follows from definitions that 

oo 

taUx- 1 } + {y- 1 }) = 

i,j=0 



di dj 



/jri 3^ 1 y Sj 



-3i—l„—3j 



X 



-1 



y 



-i 



-a(x,y) 



which by and the Fay identity 



implies that 



det 



Xi + yj _ 



{ - 1)n mM TAiUx7l]+Uy - 1}h (37) 



3.4 Asymptotics and KdV equations 
3.4.1 

Combining the formulas ( |3.2|) and (|3.7|) , we obtain the following asymptotics 



det [K {-xl-y])} 



l<i ,j<n 
1 



e,e'e{±l} n 

A(ex) A(e'y) 



+ E 



{- 



(3.8) 



n (etXi + tyj) 

where ex denotes the term-wise product of two vectors. 

Since the left-hand side in ( |3.8| ) is analytic and ( |3.8| ) is valid in the complex 
domain, it can be extended to the asymptotics of the integrand in ( |3.1|) by 
letting y { -> x { in (ggp. 

The full form of the ?/; — > limit of ( |3.8| ) is rather messy, but fortunately, 
we can concentrate only on a small fraction of the terms that survive all 



selection rules discussed in Section 3.1 



3.4.2 

Most importantly, we get nonoscillating terms only if e' = —e which cuts the 
number of summands from 4 n down to 2 n . These 2 n terms can be conve- 
niently organized using the following operators V;. 

Let f(xi, yi, . . . ) be a function of Xi, yi and some other variables which is 
supersymmetric in x; and y\ in the sense that its restriction to the diagonal 
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does not depend on x i} that is, 

= (3.9) 

Introduce the following operator 

fv7 ,-]/ x ,. l/(si,yi, •••) -/(-xi,-yi,...) 
Va/ (xa, ...) = hm . 

Note that the numerator here vanishes on the diagonal y% = x±, so this limit 
is well defined. 

With this notation, we see that 

det [K (-xl -..j) ,, r ( , ~ { _ 2ni ^ UiXi ViV 2 •••V,.<I. • ... (3.10) 
where 

• = II*.II^(E{y-E{y), (3.H) 

A = «p (§i(«» -,?)), ^ = gzgkzg, (3.12) 



and dots in ( |3.10|) stand for oscillating terms. Observe that all factors in 



( |3.11| ), ( p.!2| ) have the supersymmetry ( |3.9|) . 



3.4.3 

The operators V» commute and satisfy the following Leibnitz-like rule 

Vi(/ ■ <?) = V*(/) ■ <5 4 (s) + *(/) ■ Vi(g) . (3.13) 

Let us examine the effect of applying V, to (|3.11|) . 

Denote by r M the coefficients in the expansion of the ta 

^E^ < 314 > 

where the summation is over all partitions //, | Aut fi\ is the product of facto- 
rials of multiplicities of parts in //, = Y[ t^- The variables tj, are specialized 
in ( 3.11|) in the following way 
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Recall that ta depends only on with k odd. It is clear that 



- x ■ 



-fc-i 



k odd . 



(3.16) 



and also that Si simply removes and from the sum ( |3.15| ). 
Similarly, it is clear that 



Vi (f>i = 2ix { , 



(3.17) 



Finally, for the last type of factors in ( 3.1 1|) we have 



Si cj)ij = Sj 4>ij = 1 
which, in particular implies that 

Vi Sj 4>ij = Vj Si 4>ij = . 

It follows that the Vj's have to be applied to 4>i/s in pairs to get a nonzero 
result 



rp 2 I />-»2 

(x? - ^ 2 ) 2 ' 



(3.18) 



3.4.4 



The factor corresponding to ( 3.18 ) in the asymptotics of det [K(— Xi, — 
will be 



x- 



iC ^ I ^ 



(3.19) 



where the first term comes from the prefactor in fl3.10| ). This has a second 
order pole on Xi = Xj, which in the full asymptotics of (|3.10| ) cancels out with 



poles of oscillating terms. This singularity of ( |3.19| ) is immaterial and can 
be removed as explained in Section |3.1.4| . What is important about ( 3.19 ) 
is that is has degree —2 in Xi and Xj and, hence, whatever contribution it 
makes to the asymptotics of ( ft.ip , it will be in degree in & and Since we 
are interested in terms of strictly positive degree in all variables, see Section 
p.l.6| , the terms containing ( |3.19| ) are negligible. By the same token, any 
terms containing (|3.17|) can be also discarded. 
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It follows that relevant terms in the asymptotics ( |3.10| ) are obtained by 
applying all operators V; to the r-function. From (|3.16|) and the Leibnitz 
rule Q3.13D we obtain 



(— l) n | Aut fx\ m M , £(/i)=n, 
otherwise . 



where denotes the monomial symmetric function 



" - I Ant n\ ll^W • 

creS(n) 



Aut jJL 

Therefore, modulo irrelevant terms, the asymptotics of the integrand in ( pTT| 
is 

det [K(— X;, —Xa)] ~ r — — > 



e(fi)=n 

Recall that only partitions \x with odd parts enter this sum 
3.4.5 

Now it remains to use the formulas 

e -t* x «-l dx „ _W £_>+(), 

e 

and 

2Jfe + l\ (-2) fc+1 v^F 



r 

to obtain the expansion 



(2k + 1)!! 



mi,...,ra„ v y 

where dots stand for irrelevant terms, that is, for diagonal terms and terms of 
negative degree in the £j's, and r 2m+ i is the coefficient in ( |3.i4|) corresponding 



to the nonincreasing rearrangement of the numbers 2mj + 1. 
By ( [2.8D and Proposition |2.5| we have 

5(x)= R°(x a ) = R°(x) + ... , 
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where R° is the connected part of R and dots stand for diagonal terms. Now 
( 3.2Q) and Theorem [I] imply that 



n 



X ■ 



vmi,...,m„ 



(2mi + 1)!! 



+ . . . 



where r denotes the rescaled r-function 

^ = H) ImI 2 ImI/3 Tm) 



(3.21) 



circle stands for the connected part, and dots stand for terms of negative 
degree in x, that is, for the g = terms in the case n — 1, 2. This concludes 
the proof of the KdV equations. 

Note that the rescaling (|3.21|) is equivalent to the following rescaling 



(6&-1)!! 

36 fc (2A;)! 



keeping the relation = ^|| a k . 
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